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Effects of low-lying fermion modes in the e-regime* 
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We investigate the effects of low-lying fermion modes on the QCD partition function in the e-regime. With 
the overlap Dirac operator we calculate several tens of low-lying fermion eigenvalues on the quenched lattice. By 
partially incorporating the fermion determinant through the truncated determinant approximation, we calculate 
the partition function and other related quantities for Nf — 1 and compare them with the theoretical predictions 
obtained by Leutwyler and Smilga. 



1. Introduction 

The low energy behavior of QCD, and therefore 
lattice QCD, is effectively described by the chiral 
Lagrangian. Approaching the massless limit, the 
pion's Compton wave length becomes longer than 
the spatial extent of the lattice, and the zero mo- 
mentum modes dominate the partition function. 
This is the so-called e-regime of QCD, where var- 
ious properties of the QCD partition function is 
known analytically j^. For instance, for one- 
flavor QCD the partition function is given by 
Z = e"*^^™"^ for finite 9 at the leading order. 
Here, m is the quark mass, S is the chiral con- 
densate, and V is the space-time volume. It can 
be divided into topological sectors as 



(1) 



where Iu{x) denotes the modified Bessel function, 
and V is the topological charge. 

For lattice QCD to study the small quark mass 
region, it is therefore the first step to confirm 
such analytic results. Eventually one can deter- 
mine the low energy constants of chiral perturba- 
tion theory by measuring correlation functions in 
the e-regime, as discussed for example in In 
this exploratory study we investigate the mass de- 
pendence of the partition function, the sum rules 
for eigenvalues of the Dirac operator, as well as 
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the topological susceptibility. The effect of low- 
lying fermion eigenmodes is incorporated using 
the reweighting technique, and thus we study the 
above relations in the one-fiavor Nf = 1 case. 

Since the chiral symmetry is essential in de- 
riving the partition function in the e-regime, 
we employ the overlap Dirac operator D = 
[1 -I- 75Sgn(7Jiv)]- To compute the sign func- 
tion sgn(i?iy), 14 lowest eigenmodes of Hw are 
treated exactly and we used the Chebyshev poly- 
nomial of degree 100-200 to approximate the rest 
of the eigenmodes. On a IC* lattice at /3 = 5.85 
we calculated 50 lowest eigenvalues of D using the 
implicit restarted Arnoldi method (ARPACK). 
The topological charge is obtained from the num- 
ber of zero modes for each gauge configuration, 
and we accumulated 168, 290, and 149 configura- 
tions for \v\ = 0, 1, and 2, respectively. 



2. Truncated determinant approximation 

In order to study the effect of fermion deter- 
minant, we use an approximation to include only 
the low-lying fermion eigenmodes. Namely, the 
fermion determinant is replaced by a product 
Jlt^T^'d'^il^ -I- m^), where \i is the i-th lowest 
eigenvalue of the overlap Dirac operator D for 
a given configuration. If A'max is large enough 
to cover physically relevant eigenmodes, this ap- 
proximation is expected to reproduce the right 
physics. For our choice of parameters Aio lies 
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Figure 1. Correlation between Y[i=i l^il"^ ^^"^ 
n^=i fo'" the topological sector 0, 1, and 
2. For the non-zero topological charge, the zero 
modes are subtracted. 
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Figure 3. Partition function of a given topolog- 
ical sector divided by the sum of the partition 
function of all topological sectors. The curves are 
the analytic prediction pp. 
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Figure 2. Truncated determinant as a function of 
-^max, normalized with an average over five con- 
figurations shown in the plot (all in the = 
sector). 



around 700 MeV, and therefore should already 
be large enough to be treated as the ultraviolet 
mode. The neglected eigenmodes are the ultra- 
violet modes, which should be irrelevant to the 
low energy physics. Such idea was implemented 
previously as an algorithm for the Wilson fermion 
[3], in which the ultraviolet modes are also incor- 
porated using either the effective gauge action or 
the multi-boson method. See also 

In order to demonstrate how the truncation 
works, in Fig. ^ we plot a correlation between the 
truncated determinants with A^max = 10 and 50 



for m = 0. Each dot in the plot represents a gauge 
configuration, and the strong correlation among 
the different truncation indicates that the prod- 
uct Hi^ii l-^iP is mostly independent of gauge 
configuration. In Fig. |2] we show how the trun- 
cated determinant varies as a function of A'max- 
We can observe a rapid change below TVmax ~ 10, 
with which the order of magnitude is determined. 
The values still changes slowly below iVmax ^ 50, 
and stays almost fiat above it. 

3. Partition function 

In Fig. Owe show the mass dependence of the 
partition function Ziy{m) / Z (m)\g^Q for the topo- 
logical sectors ji'l = 0, 1, and 2 at Nf — 1. The 
truncated determinant with iVmax = 50 is used. 
We can see that the data agree with the ana- 
lytic expectation I,y{m'SV)/ exp{m'S,V) PP shown 
by the curves. From this fit we obtain E = (243 
MeV) 3. 

By differentiating the partition function with 
respect to the quark mass, Leutwyler and Smilga 
obtained a sum rule for the eigenvalues |Ij 

Left hand side of this equation is quadratically di- 
vergent in the ultraviolet region, since the eigen- 
mode distribution behaves as A"^ in the bulk. 
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Figure 4. Difference of {J2n ^/{^n^V)^)iy among 
different topological sectors. Analytical expecta- 
tion is shown by horizontal lines. 




Figure 5. Mass dependence of the topological 
charge squared (i'^) for Nf ~ and 1. 



Therefore the sum rule makes sense only af- 
ter the ultraviolet divergence is regulated and 
the particular volume dependence is extracted. 
Here, instead, we consider a difference of be- 
tween different topological sectors. Because the 
bulk distribution is expected to be independent 
of topology, the ultraviolet divergence cancels in 
such differences. Fig.0]shows these differences as 
a function of the ultraviolet cutoff. We find that 
the results are completely independent of the cut- 
off and agree with the analytical predictions. 

4. Topological susceptibility 

Topological susceptibility x is obtained from 
the partition function by differentiating with re- 
spect to 9^. On the lattice one can simply cal- 
culate the expectation value of the topological 
charge squared (z/^) = x^- In Fig- IH] we plot 
(j/^) as a function of quark mass for both Nf 
= and 1. Near the massless limit the topo- 
logical susceptibility grows linearly for Nf — 1 
as expected from the form of the partition func- 
tion e"iSVcose Pqj. larger quark masses it sat- 
urates toward the quenched value corresponding 
to X = (205 MeV)^. This behavior was previ- 
ous observed in 0. By fitting the slope near the 
massless limit we obtain E = (238 MeV)^, which 
agrees well with the determination from the par- 
tition function itself. 



5. Conclusions 

Using the truncated determinant approxima- 
tion of the overlap Dirac operator, we confirmed 
that the Leutwyler-Smilga's analytic predictions 
for the QCD partition function is satisfied for 
Nf — 1. Results are stable against the number 
of eigenmodes included, if it is greater than 10 
(Aio ~ TOOMeV). It suggests that the dynamical 
fermion algorithms designed to separate low and 
high eigenmodes, such as [H], work effectively. For 
larger volume or smaller lattice spacing the num- 
ber of eigenvalues to cover the equivalent physical 
region becomes much larger (proportional to the 
lattice volume), and their calculation could be nu- 
merically prohibitive. Still it would be helpful to 
isolate the low-lying modes, as they account for 
the dominant part of the fermion determinant. 
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